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Abstract 

We propose a model dielectric function for ropes of single-walled nanotubes dis- 
tributed in a glassy graphite host medium. We study the significance of the bosonic 
charge excitations arising in interacting quasi-one-dimensional systems in the screen- 
ing processes. We also pay special attention to the role of the intertube Coulomb 
interactions. In order to compare with experiments, weak relaxation processes are 
also considered in the relaxation-time approximation. 
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1 Introduction 



Carbon nanotubes are being considered as substitutes of silicon in many elec- 
tronic devices (see [1] for a review). Besides their extraordinary chemical sta- 
bility and strength, they support ballistic electronic transport in a large range 
of temperatures. On the basis of these properties and several experimental 
demonstrations of performance in specific devices [2], there is an ongoing ef- 
fort to find more controllable and efficient methods of synthesis, which are 
necessary in mass-production electronics. Simultaneously, the importance of 
the one-dimensional nature of the electronic structure and the role of the elec- 
tronic interactions in these materials has been recognized. The basic electronic 
properties of single-walled nanotubes (SWNT) are determined by the helicity 
of the graphene sheet that forms the nanotube, i.e. the chirality. Nanotubes 
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with (n, n) or (0, n) chiral vectors are metallic. In these systems the interac- 
tions, either Coulomb [3] or phonon-assisted [4], play a very important role 
in the transport properties of nanotubes at low temperatures. The observa- 
tion of Luttinger Liquid (LL) phenomenology [5,6] has been attributed to the 
existence of strong interactions in the forward scattering channel [7,8], along 
with the one-dimensional character of the electronic structure. In addition to 
clear signals of non-BCS superconducting behavior both in ropes [9] and in 
small-radius single nanotubes [10], other types of charge instabilities have been 
predicted theoretically [11,12,13], deriving from the strong renormalization of 
the electron-electron couplings typical of quasi-one-dimensional systems at 
very low frequencies. 

In this paper we study the screening processes involved in a rope of nanotubes 
and we compute the dielectric response. Our calculations are valid in the fre- 
quency range A < u < E c , where A is the single-particle gap associated with 
the low temperature instabilities in the rope, and E c is an energy scale, to be 
established later, below which the ropes behave as one-dimensional interact- 
ing electron systems. Here, we will be mainly concerned with the longitudinal 
response. The real part of the conductivity can be expressed in terms of the 
imaginary part of the dielectric function via Hopfield-like formulas. 



2 Dielectric response 

We start the discussion giving precise meaning to the range of the Coulomb 
interaction in nanotubes. Unlike in three dimensional metals, in strictly one- 
dimensional systems there is no conventional plasmon screening mechanism of 
the Coulomb interaction that remains long-ranged [14,15,16]. On the other 
hand, SWNT's have an extra dimension perpendicular to the main axis, 
namely the radius of the tube R, and therefore the electronic wave function has 
an extension in this direction which cuts off the otherwise diverging Fourier 
transform U (k) of the interaction. There are several ways to parametrize this 
partially screened Coulomb interaction. Here, we use one that, despite its sim- 
plicity, captures the physical mechanism of the screening at small momentum 
transfer using just two parameters, namely U(k) = 7TVF 2 Uo log (^p), and it is 
well motivated from the physical point of view [17]. Uq encodes the intensity 
of the interaction and k c is a momentum cut-off determined by a length scale 
associated to the transverse direction. 

The role of the Coulomb interaction in an isolated metallic nanotube has 
been the subject of intense study, especially in relation to the experimental 
observation of Luttinger Liquid phenomenology in these systems. Actually, it 
has been shown [7,8] that the strength of the backscattering and Umklapp 
processes mediated by the Coulomb potential is reduced by a relative factor 
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~ 0.1 (a/R), where a is the lattice spacing, with respect to forward scattering 
processes. In addition, the intertube single particle hopping is exponentially 
suppressed in typical ropes which are made of nanotubes with different diam- 
eters and chiralities [18]. Under these conditions the polarization effects can 
be computed within the RPA approximation. Let us consider then the generic 
non-interacting polarizability: 



dp f[e 8 (p)]-f[e a >(p + k)] 
2ir uj + e s (p) - e s >(p + k) 



where s and s' are arbitrary band indices, and first we estimate the inter- 
band polarization effects. We recall some elementary facts about the electronic 
band structure of (n, n) nanotubes. For nanotubes with typical diameters of 
~ 1.5 nm, curvature effects are negligible and the band structure can be ob- 
tained by wrapping the graphene sheet and finding quantization conditions 
for the wave vector in the radial direction. One finds 
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where t n ~ 2.5 eV is the hopping integral between two carbon ir orbitals, 
L y is the perimeter of the tube measured in number of hexagons, and s is a 
band index (s = 0, ±1, . . .). Therefore, there are two bands (s = 0) crossing 
the Fermi level at exactly two Dirac points q = ±k F = ±-^^. Linearizing 
the dispersion relation, Eq. (2), we have e s=0 (p) = ±vp\p ± k F \. The next 
completely empty (and full) bands s = 1 (s = — 1) have two extrema at 
p ~ ±kp. We can expand e s (p) as e s (p) ~ +E- int +p 2 /2m, where, according to 
Eq. (2), E int = ±t n ^j2 — 2 cos(n/L y ). For the first empty band E- int is of the 
order of 1.5 eV for nanotubes with typical diameters (~ 1.5 nm). For momenta 
close to ±&f the interband polarizability is: 
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and only contributes at large frequencies. 



On the other hand, there is a remarkable intertube screening effect. Let us 
call g and u the bare intratube and intertube coupling constants, respectively, 
and D and V the corresponding vertices. Since the intertube polarizability 
between two metallic nanotubes is of the same order of the intratube one, 
the effective density- density interaction between two electrons located in the 
same nanotube is corrected by the interaction of electrons in different 
nanotubes V^\a,b). We are using the conventional notation for super and 
subindices [19]. The other pair of indices, a and b, label different nanotubes. 
We consider spin independent interactions = g^2 and we define our bare 
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couplings as gf ' = g\V + gfl (analogously for the w's). The virtual processes 
that renormalize the interactions can be represented in diagrammatic form 
(see Fig. 1) and translate into the following Dyson equations for the vertices: 



d? = g? + E ^wPnoff 

klmn 

+ E Y. r lkraiTlinjU k l \a,c)I\y^\c,a) (4a) 

klmn c^a 

Vf\a, b) = u^(a, 6) + E E VkmimnjU k l \^ c)UVi n \c, b) 

klmn Cy^a,b 

+ E VkmiVm^ia, &)£>£> + gPv£>(a, b)} (4b) 

klmn 

where r) kmi = 5 4i S km + 5 2i (l - S km ); i] lnj = 5 Aj 5i n + 5 2i (l - 5 ln ) and n (0) = 
5jlU f + (i _ fy)nL°> , 5 rv are Kronecker 5 functions, and nf and nL 0) are 
the respective polarizations for right and left branches with linear dispersion, 
given by 

n f(fc,u;) = -— -; n W (A:,cu) = - . ~^ - . (5) 

+ v ; 27r(w-t; F A;)' v ; 2tt(lu + v F k) y ' 

superscript zero refers to the pure case. 

For a single nanotube we can decouple the equations in the form 
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where we used the combinations ^ ± ^ = 
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Using the 



Eqs. (5) we obtain the expressions for the vertices. For instance, in the 
positive branch is: 
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Fig. 1. (a) Diagrams contributing to the screening of Coulomb interactions between 
currents l a within the same nanotube. The currents have well-defined chirality and 
spin. The last term takes into account the polarization of the rest of n — 1 metallic 
nanotubes in the rope, (b) Diagrams contributing to the screening of Coulomb 
interactions between currents l a and lb belonging to different metallic nanotubes. 



In the limit of frequencies larger than the "plasma" and free electron energies 
(uj ^> u±k,v F k) we recover the value of the bare coupling d\^(uj — > oo, k) — 
gf \ If all the bare intratube couplings g\^ = U(k), the intratube vertex Eq. (7) 
takes the form 

from which one can extract the dielectric function as (D*f* = U/e^): 



nv F uj 1 — [v F ky 

This result can be obtained by using the alternative route of finding the ele- 
mentary bosonic excitations using the bosonization technique, and after that 
computing the density-density correlations. 

Let us now consider the problem of a bundle of n nanotubes. We assume that 
the perpendicular size of the rope is small and therefore, in first approxima- 
tion, electrons in different nanotubes interact with the same bare strength 
v!f\a,b) = u , whereas = g . Then Eqs. (4) for the vertices can be rear- 
ranged in compact matrix form as e^ ) (D V) T = (go u ) T , where the 2x2 



matrix 

e(0) = /l - 2 flb n<°> -2(n-lK#) \ 
{ -2 Mo n(°) l-2[g + (n- 2)u ]U^ ) 

generalizes the dielectric function to the case of intratube and intertube screen- 
ing, and n(°) = nf + rrl 0) . 

The two eigenvalues of e^°\ Eq. (12), give a measure of the screening due to 
both intratube and intertube scattering. We assume the same functional de- 
pendence for both intratube and intertube couplings, go = U\(k) and uq = 
U2{k), both given by U(k) but now characterized by different intensities, Uq\ 
and U02, and different cutoff momenta, k cl and k c2 , respectively. In particular, 
different cutoff momenta can be related to the different radii of a single nan- 
otube and of a nanotube rope. In that case, the eigenvalues of Eq. (12) can 
still be expressed analytically in terms of the total polarizability U.^°\ and will 
be discussed numerically in the presence of impurities and of a host medium 
below. 

The low energy sector of the theory is integrable and therefore we expect fi- 
nite weight for the Drude delta peak in the conductivity at zero frequency. 
In order to make contact with experiments, we need to introduce weak relax- 
ation effects, due to collisions with defects or impurities. We can introduce 
these effects in the non-interacting polarizability within the relaxation time 
approximation, which is equivalent to carry out the substitution [20,21]: 

n<°>r* U ) - (i + »Hn (0) (^ + »/r) 

1 ' ' l + (i/u;T)UW(k,u; + i/T)/U^(k : Oy { } 

to the total polarizability Il( ) = 11^ + U^. Here, r is an external relaxation 
time which can be extracted phenomenologically by comparison and fitting of 
experimental results [22,24]. In the case of a single nanotube, Eq. (11), or for 
a rope of n nanotubes with intertube interaction simply proportional to the 
intertube one U mt (k) = r]U(k), we explicitly find 

4 iiW) = 1 _ &ok *(^)_£_, (14) 

where momenta are measured in units of (fpr) -1 and frequencies in units of 
r _1 , and U — [1 + r)(n — l)]Uo. In the case of a rope of n nanotubes, with 
intratube and intertube interactions parametrized by U\{k) and U 2 (k), respec- 
tively, one has to perform Mermin's substitution for the total polarizability, 
Eq. (13), in the elements of the dielectric matrix, Eq. (12), before solving the 
eigenvalue equation. Fig. 2 shows the real and imaginary parts of the larger 
dielectric eigenvalue, as a function of k (left panels), and of u (right panels). 

Our results are consistent with those found in semiconductor quantum-wire 
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Fig. 2. Momentum (left panels) and frequency (right panels) dependence of the real 
(upper panels) and imaginary (lower panels) parts of the larger dielectric eigenvalue 
of Eq. (12), including weak relaxation effects. Momentum and frequency are rescaled 
as vprk and lot, respectively. We take n = 10, with J7oi = 7, J7o2 = 2, vprkd = 0.25, 
and vprk C 2 = 0.05. 

nanostructures [14,23]. Specifically, a system of two arbitrarily close wires with 
equal densities and effective masses is formally equivalent to the single nan- 
otube we have considered. In our case, = U (k) and the acoustic plasmon is 
degenerated with the continuum. However, strong renormalization effects are 
expected at low frequencies [11,12,13] and it would be interesting (but beyond 
the scope of this communication) to study the impact of such effect in the 
plasmon peaks, as observed in far infrared spectroscopy or Raman scattering. 
Likewise, small-radius nanotubes seem good candidates for the observation 
of three plasmons and an analysis of plasmon stability and spectral weight. 
Following Ref. [23], one can define the dispersion relations associated with 
the RPA collective modes as the zeroes of det e in Eq. (12), now taking into 
account also the effects of impurities. One finds two plasmon-like modes, lin- 
early dispersing in the long- wavelength limit. Specifically, the lower- frequency 
mode becomes more damped with increasing number n of nanotubes in a 
rope, which is consistent with the results of Ref. [23], where the lower (acous- 
tic) mode approached the Landau damping region with decreasing spatial 
separation between the two wires. 

So far, we have considered the intrinsic properties of a nanotube rope. In 
realistic samples the nanotubes are scattered in a medium of glassy graphite 
produced during the synthesis. Following [24] we will now assume that the 
nanotubes are placed in a host medium of insulating glassy graphite with a 




Fig. 3. Real (upper panel) and imaginary parts (lower panel) of total dielectric 
response, Eq. (15), corresponding to the larger eigenvalue of the dielectric matrix, 
Eq. (12), as a function of wr, for vprk = 0.1, = 1, and different values of the 
filling factor /. All other parameters are as in Fig. 2. 

dielectric constant e h . The dielectric reponse of such a composite is given by 
a generalized Maxwell-Garnett expression [25] 



e T 2(l-/) + (l + 2/)(e/e h ) 
e h ~ (2 + /) + (l-/)( c / Cfc ) : 



(15) 



where / is the filling fraction of the nanotubes. Expression (15) interpolates 
between the the host dielectric function at / = and the intrinsic value 
e(u), given by Eq. (14), at / = 1. The value of the host medium can 
be computed within a classical dispersion theory [24]. We assume that is 
nearly constant in the frequency range we are interested. The most valuable 
information that can be extracted when comparing with real experiments is 
the position and width of the maximum in the conductivity at low frequencies 
(which is simply related to the dielectric constant). In Fig. 3 we show the total 
dielectric response in a medium with = 1 as a function of out for different 
values of the filling fraction /. 

In conclusion, we have studied the nature of the screening processes of metallic 
single-walled nanotube ropes. We have considered different factors contribut- 
ing to the longitudinal dielectric response of these systems: intratube and in- 
tertube Coulomb interactions, the presence of a glassy graphite environment 
and the influence of weak relaxation effects produced by impurities or defects. 
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Our results suggest that metallic nanotubes may be a new playground for the 
study of collective charge excitations 
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